Spin-Hall Effect in Two-Dimensional Electron Systems with Rashba Spin-Orbit 

Coupling and Disorder 
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Using the four-terminal Landauer-Biittiker formula and Green's function approach, we calculate 
numerically the spin-Hall conductance in a two-dimensional junction system with the Rashba spin- 
orbit (SO) coupling and disorder. We find that the spin-Hall conductance can be much greater or 
smaller than the universal value e/8n, depending on the magnitude of the SO coupling, the electron 
Fermi energy and the disorder strength. The spin-Hall conductance does not vanish with increasing 
sample size for a wide range of disorder strength. Our numerical calculation reveals that a nonzero 
SO coupling can induce electron derealization for disorder strength smaller than a critical value, 
and the nonvanishing spin-Hall effect appears mainly in the metallic regime. 

PACS numbers: 72.10.-d, 72.15.Gd, 71.70.Ej, 72.15.Rn 



The emerging field of spintronics,Q,0] which is aimed 
at exquisite control over the transport of electron spins 
in solid-state systems, has attracted much recent interest. 
One central issue in the field is how to effectively gener- 
ate spin-polarized currents in paramagnetic semiconduc- 
tors. In the past several years, many works 0, 0, 0, |f| 
have been devoted to the study of injection of spin- 
polarizcd charge flows into the nonmagnetic semicon- 
ductors from ferromagnetic metals. Recent discovery of 
intrinsic spin-Hall effect in p-doped semiconductors by 
Murakami et al. and in Rashba spin-orbit (SO) cou- 
pled two-dimensional electron system (2DES) by Sinova 
et al. 7] may possibly lead to a new solution to the is- 
sue. For the Rashba SO coupling model, the spin-Hall 
conductivity is found to have a universal value e/87r in a 
clean bulk sample when the two Rashba bands are both 
occupied, being insensitive to the SO coupling strength 
and electron Fermi energy 0. 

While the spin-Hall effect has generated much inter- 
est in the research community, [111 ES El 13 13 El 

Il5l Hil Kl\ theoretical works remain highly controversial 
regarding its fate in the presence of disorder. Within a 
semiclassical treatment of disorder scattering, Burkov et 
al. E3 and Schliemann and Loss E2 showed that spin- 
Hall effect only survives at weak disorder. On the other 
hand, Inoue et al. E3 pointed out that the spin-Hall ef- 
fect vanishes even for weak disorder taking into account 
the vertex corrections. Mishchenko et al. EH further 
showed that the dc spin-Hall current vanishes in an im- 
pure bulk sample, but may exist near the boundary of a 
finite system. Nomura et al. evaluated the Kubo for- 
mula by calculating the single-particle eigenstates in mo- 
mentum space with finite momentum cutoff, and found 
that the spin-Hall effect does not decrease with sample 
size at rather weak disorder. Therefore, further inves- 
tigations of disorder effect in the SO coupled 2DES are 
highly desirable. 



In this Letter, the spin-Hall conductance (SHC) in a 
2DES junction with the Rashba SO coupling is studied by 
using the four-terminal Landauer-Biittiker (LB) formula 
with the aid of the Green's functions. We find that the 
SHC does not take the universal value, and it depends 
critically on the magnitude of the SO coupling, the elec- 
tron Fermi energy, and the disorder strength. For a wide 
range disorder strength, we show that the SHC does not 
decrease with sample size and extrapolates to nonzero 
values in the limit of large system. The numerical cal- 
culation of electron localization length based upon the 
transfer matrix method also reveals that the Rashba SO 
coupling can induce a metallic phase, and the spin-Hall 
effect is mainly confined in the metallic regime. The ori- 
gin of the nonuniversal SHC in the 2DES junction is also 
discussed. 

Let us consider a two-dimensional junction consisting 
of an impure square sample of side L connected with 
four ideal leads, as illustrated in the inset of Fig. 1. The 
leads are connected to four electron reservoirs at chem- 
ical potentials fix, fi2 and ^3. In the tight-binding 
representation, the Hamiltonian for the system including 
the sample and the leads can be written as [lil 

H = -t C la C 3,<r +^£ict ff C iCT 
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Here, V so is the SO coupling strength, = in the leads 
and are uniformly distributed between [—W/2,W/2] in 
the sample, which accounts for nonmagnetic disorder. 
The lattice constant is taken to be unity, and S x and 
8 y are unit vectors along the x and y directions. In the 
vertical leads 2 and 3, V so is assumed to be zero in order 
to avoid spin-flip effect, so that a probability-conserved 
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spin current can be detected in the leads. 

The electrical current outgoing through lead I 
can be calculated from the LB formula 20] Ii = 

(e 2 /M Ejvi T l,l'( U l' - U 0, where U t = m/(-e) and T u , 
is the total electron transmission coefficient from lead I' 
to lead I. A number of symmetry relations for the trans- 
mission coefficients result from the time-reversal and in- 
version invariance of the system after average of disorder 
configurations, use of which will be implied. We consider 
that a current I is driven through leads and 1, and ad- 
just Ui's to make h = — Jo = / and h = h = 0. Since 
in the present system the off-diagonal conductance G xy 
vanishes by symmetry, Uo — U\ equals to the longitudinal 
voltage drop caused by the current flow /. In the vertical 
leads 2 and 3, where V^ = 0, the electrical currents are 
separable for the two spin subbands // = J/f + ijjj with j 
and I for spins parallel and antiparallel to the z-axis. The 
spin current is given by 1^ — [K/2(—e)] (/;•[■ — IjjJ. By 
use of the LB formula, it is straightforward to obtain for 

(3) (2) 

the transverse spin current I sI ^ = —I^g = G s h\Uq — Ux). 
Here, the proportional coefficient 



G 



sH 



(2) 



is the SHC, where T3er,o is the electron transmis- 
sion coefficient from lead to spin-cr subband in lead 
3. Equation can be calculated in terms of the 
nonequilibrium Green's functions [2ll E2I G s h — 
-(e/47r)Tr(r 3 77G"T G a ). Here, 77 = 1 and -1 in the 
spin-| and spin-J, subspaces, respectively, and Ti = 
— (£;)'] with S; the retarded electron self-energy 
in the sample due to electron hopping coupling with lead 
I. The retarded Green's function G r is given by 



G r = 



1 



E-H a - Ef= (Si) 



(3) 



and G a = (G r )t, where E stands for the electron Fermi 
energy, and H c is the single-particle Hamiltonian of the 
central square sample only. The self-energies can be first 
computed exactly by matching up boundary conditions 
for the Green's function at the interfaces by using the 
transfer matrices of the leads [24|. The Green's function 
Eq. J2J is then obtained through matrix inversion. In 
our calculations, G s h is always averaged over up to 5000 
disorder realizations, whenever W =/= 0. 

In Fig. 1, the SHC G s h is plotted as a function of the 
electron Fermi energy E at fixed size L = 40 for several 
disorder strengths. The SHC is always an odd function 
of electron Fermi energy E, and vanishes at the band 
center E = 0. The antisymmetric energy dependence 
of the SHC is similar to that of the Hall conductance 
in a tight-binding model (2^|, and originates from the 
particle-hole symmetry of the system. For E < and 
E > the charge carriers are electron-like and hole-like, 
respectively, and so make opposite contributions to the 
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FIG. 1: Spin-Hall conductance G s h for some disorder 
strengths as a function of electron Fermi energy E. Here, the 
sample size L — 40 and the spin-orbit coupling Vso ~ 0.5t. 
Inset is a schematic view of the four-terminal junction. 
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FIG. 2: Spin-Hall conductance as a function of sample size L 
for different disorder strengths at E = — 2t and Vso = 0.5t. 
Error bars due to statistical fluctuations, being smaller than 
the symbol size, are drawn inside the open symbols. Inset: 
spin-Hall conductance as a function of disorder strength for 
L = 40 and 60. 



SHC. With increasing E from the band bottom E ~ — At, 
except for a small oscillation due to the discrete energy 
levels in the finite-size sample, G s h increases continu- 
ously until E is very close to the band center E = 0. It is 
easy to see from Fig. 1 that at weak disorder W % t the 
calculated G s h may be greater than the universal value, 
namely, 0.5 in our unit e/4ir. 

In order to determine the behavior of the spin-Hall ef- 
fect in large systems, we calculate the SHC as a function 
of the sample size from L = 10 up to 100 for different 
strengths of disorder, as shown in Fig. 2. For weak dis- 
order W < 3i, the SHC first increases with increasing 
sample size, and then tends to saturate. In particular, 
for W < t, we see that the SHC can be several times 
greater than the universal value e/87r, when the system 
becomes large. For a stronger disorder 3i < W < 5t, the 
SHC is roughly independent of the sample size, and ex- 



3 




FIG. 3: Spin-Hall conductance as a function of spin-orbit cou- 
pling strength for some disorder strengths. Here, the sample 
size L = 40 and the electron Fermi energy E = —2t. 



trapolates to a finite value in the large-size limit. There- 
fore, it is evident that the SHC will not vanish in large 
systems in the presence of moderately strong disorder 
W < 5i. With further increase of W, the SHC becomes 
vanishingly small at W > 6t, as seen more clearly from 
the inset of Fig. 2, indicating that very strong disorder 
scattering would eventually destroy the spin-Hall effect. 

We further examine the dependence of the SHC on the 
strength of the SO coupling. As shown in Fig. 3, over- 
all, the SHC increases with increasing V so in the range 
< V so < t. For W = or weak disorder, the SHC 
displays an interesting oscillation effect with a period 
much greater than the average level spacing. According 
to Eq. (J2J, the oscillation of the SHC is a manifestation 
of the oscillation of the sideway spin-resolved transmis- 
sion coefficients. For a two-terminal junction with the 
SO coupling, similar oscillation with finite sample size 
has previously been observed for the spin-resolved trans- 
mission coefficients |l9j , where the oscillation period was 
discussed to be the spin precession length L sp . If we ap- 
ply the same condition L = nL sp with n an integer and 
notice L sp ~ ivt/V so , |l9l | we can obtain for the equiv- 



alent period in the SO coupling SV so ~ nt/L. For the 
parameters used in Fig. 3, SV so ~ 0.08i, which is very 
close to the period as seen in the figure. This indicates 
that the oscillation of the SHC is due to a spin preces- 
sional effect in finite-size systems. Experimentally, V so 
can be varied over a wide range by tuning a gate volt- 
age [2(1, Ull j an d so this oscillation effect may possibly be 
observed directly. 

Electron derealization is a crucial issue for under- 
standing electron transport properties in the 2DES, and 
has already been studied experimentally by use of mag- 
netoresistance measurements j^. For this reason, we 
investigate numerically whether the Rashba SO coupling- 
can induce a universal electron derealization in the pres- 
ence of disorder. According to the well-established trans- 
fer matrix approach, [islEi^ we calculate the electron lo- 
calization length £ on a bar of essentially infinite length 
(5 x 10 5 ) and finite width L. In Fig. 4a, the normalized 
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FIG. 4: Normalized localization length as a function of dis- 
order strength calculated on long bars of length 5 x 10 5 and 
widths L = 8, 16, 32 and 64. 



localization length £/L is plotted as a function of disor- 
der strength for V so = 0.5i and L = 8, 16, 32 and 64. At 
weak disorder, £/L increases with L, indicating that the 
localization length £ will diverge as L — > 00, correspond- 
ing to an electron delocalized metallic phase. With the 
increase of W, £/L goes down and all the curves cross 
at a point (fixed point) W = W c ~ 6.3i, where £/L be- 
comes independent of bar width L. For W > W c , £/L 
decreases with L, indicating that £ will converge to finite 
values as L — > 00, corresponding to an electron localized 
insulator phase. Thus the fixed point W = W c is the 
critical disorder strength for the metal-insulator transi- 
tion. Our result is consistent with the earlier calculation 
by Ando 0]> where a metallic phase was established at 
the band center E — for a strong Rashba SO coupling. 
Here, we also study weak SO coupling. In Fig. 4b, we 
plot the result for a SO coupling strength much smaller 
than the electron hopping integral, i.e., V so = O.lt, and 
similar phase transition is also revealed at W c ~ 4.6£. 
In general, we have performed calculations in the whole 
range from strong to weak SO coupling (details will be 
presented elsewhere), and found that electron dereal- 
ization occurs for any nonzero SO coupling strength as 
the magnitude of the disorder varies. Our result is in 
agreement with the perturbative calculation of weak lo- 
calization. 30] As y so reduces, the critical W c decreases, 
and the size- independent critical £/L increases (so does 
the critical longitudinal conductance G xx [2^, H3). In 
the limit V so — > 0, we have W c — > and all electron 
states become localized, recovering the known regime of 
the two-dimensional Anderson model for electron local- 
ization j2^| . The fact that the critical £/L changes with 
V so indicates that the SO coupled 2DES belongs to the 
universality class of two-parameter scaling |3l| . Compar- 
ing W c = 6.3i calculated in Fig. 4a for V so = O.bt and 
E = — 2t with the SHC shown in Fig. 2 for the same pa- 
rameters, we see that nonvanishing spin-Hall effect exists 
mainly in the metallic regime. 
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Our numerical study addresses the spin-Hall effect in 
a finite-size junction system with leads. A comparison 
between the spin-Hall effect and the quantum Hall effect 
(QHE) can shed some light on the nonuniversal SHC ob- 
tained. For a QHE system, delocalized states exist at 
the centers of the discrete Landau levels, which are sep- 
arated by mobility gaps consisting of localized states. In 
the unit of conductance quantum e 2 /h, the Hall conduc- 
tance is known to be a sum of the topological Chern 
numbers of all the occupied delocalized states below the 
Fermi energy |25|. If the Fermi energy lies in a mobility 
gap, the Hall conductance is well quantized to an inte- 
ger. If the Fermi energy is at a critical point, where 
a delocalized state exists, the Hall conductance intrin- 
sically fluctuates between two integers. Similarly, the 
SHC is also related to corresponding topological num- 
bers of the occupied delocalized states. However, in the 
present spin-Hall systems, the delocalized states consti- 
tute a continuous spectrum without mobility gaps (or en- 
ergy gaps ^3)- -Due to the lack of a mobility gap around 
the Fermi energy, the SHC can fluctuate and does not 
show quantized plateaus. As a matter of fact, the uni- 
versal value e/8ir predicted for clean bulk systems Q is 
0.5 instead of an integer in the unit of spin conductance 
quantum e/4n (here the electron charge e in the conduc- 
tance quantum e 2 /h needs be replaced with electron spin 
h/2). For the above reason, one could not expect that 
the SHC will not change to different values under differ- 
ent boundary conditions. In the present junction system, 
the open boundary, i.e., the connection of the finite-size 
sample with the much larger semi-infinite leads is quite 
different from the essentiall y cl ose boundary used in pre- 
vious calculations 0, 0, E3, 113, E3, E3, IlE , which 
is likely the cause for the SHC to be possibly greater 
or smaller than e/8n depending on the electron Fermi 
energy, the disorder strength and the magnitude of the 
SO coupling. Notably, the analytical calculation also 
indicates that the contacts between a sample and leads 
could enhance the generation of spin currents. Our calcu- 
lations provide an important evidence that the proposed 
intrinsic spin-Hall effect 0,0 may be realized experimen- 
tally in junction systems in the presence of disorder. 

Note added: After initial submission of this paper, 
we became aware of a couple preprints by Nikolic, Zarbo 
and Souma and by Hankiewicz et al. [32j, where similar 
LB formula calculations were carried out. Despite dif- 
ferent parameter values used, their results of nonuniver- 
sal SHC robust against disorder scattering are consistent 
with ours. 
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